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On Mildly Nonlinear Partial Difference Equations of 

Elliptic Type 1 

Lipman Bers 2 

The use of the finite differences method is in solving the boundary value problem of 
the first kind for the nonlinear elliptic equation A<£= F (x,y } <j>,<j> x ,<f> v ) is justified by first show- 
ing that the problem of the corresponding difference equation has a unique solution, and 
then that the solution of the difference equation tends to that of the differential equation 
when the net unit tends to zero. Also a numerical method of the Liebmann type for the 
computation of the solution of the difference equation is developed, and these results are ex- 
tended to more general nonlinear elliptic equations. 



1. Introduction 

The method of finite difference is frequently used 
for the numerical treatment of nonlinear partial 
differential equations of elliptic type. A theoretical 
justification of this method, however, seems to exist 
in the literature only for the case of linear equations. 
In this note such a justification is given for the 
simplest nonlinear elliptic equation 



A4>=F(x,y,<l>,<l> x ,<l>y). 



(1) 



We assume that the partial derivative F^, is non- 
negative and the derivatives F$ x , F^ uniformly 
bounded, and approximate the differential equation 
by the difference equation 

<t>(x+h y y) + <t>(x,y+h) +<f>(x-h,y) +4>(x,y—h) — 4<ft(s,y) 

h 2 



'-f\v,v, 



4>{x,y) y 



4>(x+h,y)—4>(x—h t y) 
2h 

<t>(x,y+h) — (l>(x,y—h) ) 
2h j 



(2) 



We shall show that the first boundary value problem 
for this difference equation possesses a unique solu- 
tion that can be computed by a Liebmann iteration 
method, and we shall estimate the difference between 
this solution and the solution of the boundary value 
problem for the differential equation. 

These results extend almost at once to the case 
when a more sophisticated difference equation is 
used to approximate (1), and also, due to a recent 
result by Motzkin and Wasow, to the elliptic equa- 
tion 

a(x,y)<l) xx +2b(x ) y)(l) X y+c(x ) y)(l) y y=F(x J y f <l> ) (l) X) <l> y ). 

Extensions to higher dimensions and to other nets 
are also immediate. On the other hand, our method 
is tied to the maximum-principle; it does not apply 



1 The preparation of this paper was sponsored (in part) by the Office of Naval 
Research. 

2 Present address New York University, New York, N. Y. 



to such important quasi-linear equations as the 
equation of minimal surfaces or the equations of gas 
dynamics. 

We state explicitly that the reasoning reproduced 
below involves only a straightforward application of 
familiar arguments. In section 3, in particular, we 
paraphrase well-known facts in a form suitable for 
our purposes. 

2. Notations 

We denote by Q a fixed bounded domain in the 
(x,y) -plane. is the closure and 12' the boundary of 
Q. F(x,y,z,p,q) denotes a fixed continuous function 
defined for (x,y)ett and all z, p, q. We assume that 
the partial derivatives F z , F P , F q exist and satisfy 
the inequalities 

F 2 >0, (3) 



Set 



\F,Wt\<A<+~. 
Lfo^At-FfayM,,*,,), 



(4) 
(5) 



where A is the Laplace operator {&4>=4>xx J c4>vv)- Let 
v(x,y) be a fixed continuous function defined on 12. 
Boundary value problem P consists in finding a func- 
tion <l>(x,y) continuous on and twice continuously 
differentiate in ft such that 



L[<j)] = in ft, (j)=v on ft'. 



(6) 



It is known that this problem has at most one 
solution, and that a solution exists under appropriate 
smoothness hypotheses. 

We approximate L by a difference operator L h 
defined (for A>0) by 



L*[*] = A»[0]-F{a;,y,0,Z?* i ,[0] > Z?^ lf [0] }, (7) 



where 



&n[<t>(x,y)]= 

h 2 

(8) 
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D»,x[<l>(x,y)]-- 



4>{x+h,y) — 4>{x—h,y) 
2h 



z\J»(»,y)H * ( '' y+&) ^ (a>y " &) - 



(9) 



holds jor every junction <j> defined on 12 A . 
We observe first that the equation 



;/ = l 



Let P„ be the point with the coordinates (x ,y ). ma y be written in the form 
The points (x +h,y ), (x ,y +h), (x —h,y ), (x a ,y —h) 
will be called the A-neighbors of P and will be de- 
noted by Poi, • ■ •, -Po4- A- lattice domain &„ is a <£.= 
set of points, P lt . . ., P N , situated in ft, having 
coordinates that are integral multiples of h, and such 
that all points P iv , i=l, . . ., N; v=\, . . ., 4 where 
belong to 0. Neighbors of points of il h that are not 
themselves points of Q h form the boundary Q' k of tt h . 
We denote the points of 0^_by Pjt+i, • ■ • > -FVj 
and the union of ft„ and 0* by Q h . If ft is fixed, and 
X is a function defined on ft , we denote the value of \ 
at a point P 4 (or at a point P iv ) by xt (or x*)- We Since 
also set 

L h [4> t }=^{<t>t]-F{xi,yi,<t>u"D>,A<t>i],Dh,v[^\}, (7') 
where 



i=l, 



iV, 



(17) 



(18) 



x<4(i + §««) 

x--J(i^) 



X">^(l-»)>0, 



xS(>+>*> 

x-j(>-J*> 

Sx"=i, 



A»[«,]=pjZj *»-4*«|» (8') 



D/,,x[<l>i]= 'fll "' Mi. »[#«] = 



(18) implies that 

<max (0 fl , . . ., <t>n) if 4><>0 
;min (<j> a , . . ., <t> t4 ) if 4> t <0 



(- 1 



(19) 



2A 



2ft 



(9') 



The boundary-value problem P h (for some fixed 
12*) consists in determining a function <f> defined on 
12* such that 



£»[*,] = 0, i=l, . . . ,N; 

<t>j=vj, j=N+l, . 



Af. (10) 



From now on we assume that h is so small that at 
least one 12* exists, and that 

Ah<2B, O<0<1, (11) 

where is some fixed number. 

3. An Inequality 

In this section we recall some properties of the 
linear difference operator 

In^i^Anl^ + aiDn^^ + ^iDn^^-y^ (12) 



with 



7i>0. 



(13) 
(14) 



and the equality sign in (19) holds only if <t>u= . . 
= <t>U) 7f<£f = 0. 

Let m and M be the minimum and maximum of a 
solution of (17). If M>0 and <t>=M at a point P, 
of 12j>, then ^>=Mat all neighbors of P iy at all neigh- 
bors of neighbors, etc. Hence <j>=M at a point of 
12*. A similar argument holds if ra<0. Thus 

min(O,min n ^)<0i<max(O,maXi2^), i=l, . . .,N, 

(20) 
and, in particular, 

maxoj0|<maxo;M. (21) 

Consider now the boundary value problem 

h[(t>i]=ki, i=l, . . ., N; 

**=** i=iST+l, .. ., M. (22) 

It is equivalent to the system 

N M 



8 = 1 j=N+l 



i=l,...,N, 



In particular, we want to establish the existence of where the matrices (L is ) and (a ij ) depend only on 

a constant l h and 12/,. In view of (21) the homogeneous system 

C=C(A,6,d), (15) (fc = 0, v=0) has only the trivial solution. Hence 

the "Green matrix" (G ij ) = (L i /)~* of l h in 12* exists, 
depending only on A, 6 and the diameter d of 12 such and the unique solution of (22) is given by 
that the inequality 



maxn A | <t> | < C maxo A | l h [4>] | + max Q - 1 <j> | (16) 



N M 

8=1 j=N+l 
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where and since for i=l, . . ., N 

N 



N 



we call (T iJ ) the "Green boundary matrix" of l h . ( x._T. =: y-> QiSy'l >q 

Thus every function defined on il h satisfies the *=i 

identity ~ ~ 

Thus G il <G il , or \G il \<\G il \ so that 

*i=gfi w, &[*.]+ jp^Vi, t=i, . . .,iv, Ci <^ (28) 

and hence the inequality Let $(x,y) be a continuous positive function defined 

max O4 |*|<c 1 maxo Jfc |Z Jfc [0]|+c 2 maxo i |*| (23) on for which T^]>0, and 

with niax« $ 



AT M 

d=max 2)1^*1, c 2 =max 2 |r"|. (24) 

s=1 i=w We have that 

Let be a solution of (22) with /t< = 0, Vj=5ji 



min^ h [<*>] 



-<C(A,6,d)- (29) 



(Kronecker 5). Thenar 1 By property (20) of <h = y^ 71 is Ucb 14- Y^ r^ i= 1 AT 

solutions of (17), 0<*,<1, i=*l,2, . . ., N. Thus f & Al ,J /-Wi *' ' ' ' ' ' ' 
r">0 for all £,Z. Next let <£ be a solution of (22) 

with &i = 0, v N +i= . . . ssf^ssl. Then and as the term on the left-hand side is positive, 



M M 



AT" Af 

*.= S r"== E lr«|, s |3«|ld*J<. s ?«*,. 



and by property (20) of solutions of (17), 0<*,<1. Noting that ^< 1; we obtain the inequa ii ty 
nence, 

c 2 <l. (25) N 

{mino l h [$\ } 2] l#* s | <maxj7 3>, 

Now let 4> be the solution of (22) with ki=du, Vj = 0, «=i 

so that <j>i—G il . Let tl h denote the lattice domain , . . . .. . 

obtained from Q h by removing P,. Assume that whlch im P lies tnat 

<t>i>0. In Q h we have that /»[*]=0, so that by (20) Ci<G. (30) 

6<</) f <0i, i^Z. In particular, 0<<t>i v <<t>i, *>= 1,2,3,4. 

But /*[<£,] = 1, so that From (23), (25), (28) and (30) the assertion (16) 

*! d+v \ - 4 x » f ° ll0WS ' 

4 \ 4 V ' frl z " We conclude the proof by exhibiting a function 

. &(x,y) having the desired property (without aiming 

which is absurd. Hence </>,<0, and again by (20) at the best possible value of C). Without loss of 

<£*<0fori^/. Hence generality we assume that 12 is located within the 

strip 0<x<d. 

G il <0- (26) If Ad<2, we set *(z,y) = [x- (d/2)] 2 . Then 



Now set 

7, [fa] = A»fo<] + cLiD Kx [4>t\ + PtD h , yM (27) 



A fc [*] = 2, D htX [$] = 2x-d, Z?*, f [$]=0, 



so that in £l:l h [$\>2—Ad, and of course $<d 2 /4. 

i i i ^i %i ~ - l j c j i. x ^i_ 7 Hence in this case (29) holds with C=d 2 /4[2—Ad]. 

and let 6? , P , c^, c 2 be defined as before, with l h jf Jk/>2 set 

replaced by l h . Let <£ have the same meaning as 

before, and let be the solution of the boundary m=~ tanh" 1 (* \ $(x ) y) = e* x . 

value problem # V 2 / 

C[&] =S ih i= 1, . . . , N', Then max ^<^ an d 

*-0, j-tf+l,...,*. A^]=4 e -[^^] 2 , 

&=f»<0, JU$] = 2e"* [ Smh ^ /2) ] cosh ( M fe/2), Z\ ff [$] = 0, 
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Then 



so that 



7 > 2e M* sinh (fih) Ttanh (fih/2) 



\xh 



h 



\} 



Because (tanh £)/£ is a decreasing function of J, the 
expression in the brackets is, for 0<h<29/A, not 
less than 



tanh Qxfl/A) -4.4, m 



so that 



T,M>^W 



'(^> 



Thus (29) holds again with 



<7= 



2g exp{f tanh-^^j 



^(l-fltanh- 1 ^^) 

4. Existence and Uniqueness of the 
Solution of the Difference Equation 

The inequality 
maxn^ 1 -^ 11 ! < C maxnJL ft [<*> T ] -L h [<t> n ] \ 

+m&xa' k \4> I —<t> u \ 



(31) 



holds for any two functions defined on to h . 

In fact, by the mean-value theorem, the function 
<j>=<j> J —4> u satisfies the difference equation 

lA<t>]=L h W]-L h W l ], 

where l h is defined by (12) with 
ai=—F p (Xi } yi,Zi,Pi,Zi), Si=—F q (x i) y if z i) p i ^i) J 

y i =F g (x h y i ,z i ,p i ,q i ), Zi=T t <p}+(l — t,)^ 1 , 



Vt 



2*=^ {r,(^ 2 -^ 4 ) + (l-rO(^-^)|; 



Jor every function defined on ti h . 

In particular, a solution <j> of P h satisfies the in- 
equality 



maxQ A |0|<iiL = (7#2+niaxoJfl|. 



(34) 



Now let (G is ) and (T ij ) be the Green matrix and 
the Green boundary matrix of the operator A h for 
the domain ti hj and define the continuous transforma- 



tion <t>*=T t (<j>) in the iV-dimensional (<j) h 
space by the equations 



-,<Pn)- 



N 

4>f=4>i-tJ2 l G is F{x u y u <i> u D h ^4> t ] J D hty [<j> l ]} 

s = l 

M 

- S ^v t , t=l, . . , AT, 



(35) 



where £ is a real parameter, 0<£<1. A point 

0=(<£i, . . .,<piv) taken by this transformation into 
the origin is a solution of the boundary value problem 

^ h [<l>i]—tF{x t) y i) <l> ij D htX [<l> i ] 1 D hty [<l) i ]}=O i 
i=l, . . ., N, $j=vj, j=N+l, . . ., M. (36) 

Let 8 denote the domain 0?+ . . . +<t> 2 N <NK 2 +\ 1 
By virtue of (34) no point of the boundary S' of S 
is taken into the origin by T t . Hence the degree of 
the mapping T t of S at the origin (that is the iV-di- 
mensional solid angle under which the image T t (S f ) 
is seen from the origin, divided by the measure of the 
iV-dimensional unit sphere) is well defined. Since 
T t depends continuously on t, this degree is in- 
dependent of t. For t=0 the mapping T is a transla- 

tion which takes the point c/>°, with <$= ]T) T ij Vj of 

j=N+l 

S into the origin. Hence the degree is 1. Thus for 

— * — — * 

t=l at least one point <f> of S satisfies Ti(<£)— 0, and 
problem P h has a solution. 

A somewhat longer but constructive existence 
proof will be given later. 

5. Convergence of the Solution of the Differ- 
ence Equation to that of the Differential 
Equation 

We assume now that problem P has a solution 

<$>(x,y). 

Let co (7?) be the modulus of continuity of 4>(x,y) 



the it being numbers such that 0<Ti<l. Since 

conditions (13), (14) are satisfied, inequality (16) is in 12, that is, 

applicable and yields (31). 

It follows from (31) that problem P h has at most one w(ry)=l.u.b.|^(a;,l/) — 4>{%' ,y')\ for 

solution. ~ 

Set ; (a? l »),(x / ,y / )€Q l ix-xy+{y- y y<rf. 

C 2 =m^\F(x,yfififi)\. (32) gy w pj^ m p we denote the maximum f the moduli 

Then \L h (0)\<C 2 , so that applying (31) to an arbi- of continuity and of the absolute values, respectively, 

trary function <t>=<t> J and to the function <£ n = 0, we of the second derivatives of <t>(x,y) in the closed 

obtain subdomain of Q, consisting of all points in 12 having 

a distance not less than 5>0 from 12'. mf } and 

maxQ A |0|<(7(72 + Cmaxfi fe |L^[0]|+max^|0| (33) (for a fixed ?7>0)co ( 5 2) (7]) are nonincreasing functions 
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of 8. Set 



-m v 



a->o 



Km 

8-*0 



, ,(2) 



to) 



= /.,«> 



to). 



If the second derivatives of <$>(x,y) are continuous in 
12, then m (2) <+ oo and <af (ri)-*Q for r?-*0. If (j>(x,y) 
possesses continuous partial derivatives of the third 
or fourth order, we define mf , m (3) , mf\ m (4) in a 
similar manner. Finally, we denote by co*(r?) the 
modulus of continuity of v(x,y) in 12. 

Let Q h be a lattice domain such that (i) every 
point of 12^ has a_distance not exceeding e from some 
point of 12', (ii) 12^ is contained in a closed subdomain 
of 12 having a distance not less than 5>0 from 12'. 
We denote by cj) (h ' e > 8) the solution of problem P h for 
this domain. 

Let Pj be a point of 12^. By virtue of (i) there 
exists a point (x,y) of 12' such that 

\Vj-v(x,y)\<a)*(e), \<j>j-<t>(x,y)\<u(e). 

Since <£/*•«» 8) =fy, (f>(x,y)=v(x,y), we have that 

maxoJ</>^^ 5) -0|<co(e)+w*(e). (37) 

From Taylor's formula and condition (ii) it follows 
easily that for every point P t of 12^, 

IDnA^-i^iUD^yl^-i^^mfK (38) 
so that 

\F{zi,yt,<l>i,D htX [<l>t],D h j,[4> i ] } 

- { F(x,y,4>,4>*,4>v) } t\ <2Am?h. 
Also 

|A Jk [0j-(A*) < |<2 w ? ) (A). (39) 

Since L[<£] = 0, 

rrmx Qh \L h [<t>]\ <2Amf ) h+2o>f ) (h). 

Inequality (31) now yields the estimate 



maxoj* 



(h, e, 5) . 



-<f>\<2C[mfh 

+ „P(*)] + «(f)+«*(€). 



(40) 



This relation shows that 4> can be approximated by 
^(»,«,8) ^^ an y desired degree oj accuracy. In fact, 
let r7<0 be given. We determine an e<0 so that 
co(e)+co*(e) <rj/2 and choose a closed subdomain 
12* of 12 such that every boundary point of 12* has a 
distance not exceeding e from some point of 12'. For 
every A>0 let 12/, consist of all points in the interior 
12* of 12* whose coordinates are integral multiples of 
h and whose /i-neighbors belong to 12*. Let 8 be the 
(positive) distance of 12* to 12'. We choose h so 
small that 2C[m? ) h + o>? ) (h)] < v /2. Then 



If the partial derivatives of <t>(;x,y) of the third 
order exist and are continuous in 12, we may replace 

(38), (39) by 

|Z?*.*[* J — (**) «l,l^ ilf [0j — (0 F ) *| < *mi« A 2 , (380 
|A,[0j-(A0),|<|m ( f^ (390 

and instead of (40) we obtain the estimate 

max o j0<* *- 5) -ct>\< \Cm™ (Ah 2 +2h) +co(e) +co*( € ). 

(400 

If 4>(x,y) possesses continuous derivatives of order 
four, (390 may be replaced by 

\^[<t> l ]-{^) i \<lm^h\ (39") 

so that instead of (400 wc obtain 
max fi J^^^ ) -0|<JC'(2m 5 (3) +m, (4) )/i 2 + co(e)+ w *(6). 

(40") 

If the derivatives of <t>(x,y) (of order 2, 3, 4, respec- 
tively) are continuous in 12, condition (ii) is superflu- 
ous. We may drop the superscript and subscript 8 
in (40), (400, (40") and obtain, for h->0, c->0, 
uniform convergence of <t> (h,e) to <£. In fact if all 
points of Q' h are on 12' the terms with e disappear 
and we have, say in the case m (4) <+°°, that 
|* W -*|=0(A J ). 

All preceding estimates involve (a), some a priori 
estimates for the solution of P, (/3) the exact solution 
of P h . Without some a priori information on the 
solution of the differential equation, however, its 
approximation by the solution of the difference equa- 
tion can never be estimated. On the other hand, 
the difference between an exact and an approximate 
solution of the difference equation (and only approxi- 
mate solutions can be obtained in practice) can be 
estimated at once by means of (31). 

6. Solution of the Difference Equation 
P by Iterations 

In this section we describe an effective method for 
solving problem P h , and at the same time obtain a 
constructive existence proof. 

Define (for i= 1,2, . . ., N) 



/^;r 1 ,f 2 ,f 3 ,f 4 )=^+^(^,2/i, 



f 1 - 



2h 



f 3 r 5 

; — 



2h 



')• 



Then 






is nonnegative, so that there exists functions 



(41) 

(42) 

(43) 
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such that 



ytto,(z;r>,...,rtjfs...,f*i-z. 



(44) 



Also set 



««(rs • • -,f 4 )=^i (Jijr;!- 1 , • • .,r 4 } (45) 

Let <j> be a function denned on 0/, and satisfying 
the boundary condition fa=vj, j=N-\-l, . . ., M. 
It satisfies the equation L h [<j> t ]=0, i=l, . . . , N 
and is a solution of P h if, and only if, 



tives are evaluated at a point (f 1 , . . ., £ 4 ), 

r-r^^+a-Ti)??;-", 0<r,<l. 

By (41), (44), (45), 






= 1, 






i 4 

ji(<t>i',<t>iU ' - 'j<t>ii) = iSfe 
4 „=i 

that is, if 



;=i,...,iv, 



i=l,.. .,iV. (46) 



df 1 df 3 8 " bf 2 df 4 8 " 

df 4c>z i ~dr' 



so that by (42) 



d*« 4V 2*7 



Now let # (0) by any function in 14 satisfying the 
boundary condition and let the functions <£ (1) , 
<£ (2) , . . ., be defined by either of the two following 
iteration schemes. (A), set 



w 



14> 2 F 



W 



l+Z*. 



<*><" 



■s t (<t>^, . . .,^<j>) f 



i=l, 



,iV (47) 



*«»+» = #•> = 



:»,, 



j = N+l, 



M. 



(B), set 

(»+i) =Sl (^n> > _ _ -; ^(„, ); i = i t . . ., N, (48) 

where 

_ Jfc+*> if P t ,=Pj and j<i, 



ds _ 

Hence 
and setting 



. i( i+ H * 10+H 



!+?>. 



df* 






i+¥*. 



i+ 



K^' 



t^/; } otherwise, 



=* if i=iV+l, . . .,M. 
If the sequence { <£ (n) } converges, then the limit e m ^ 



it follows from (49) that 



&r 






m w+1 =max fi J(/> (w+1) -(A (w) | 



function is a solution of P h , for s f is a continuous 
function of its arguments. 

In the case of the Laplace difference equation, (A) we conclude from (50) that 
corresponds to the difference equation of heat con- 
duction whose solution converges to that of the 
Laplace equation for n— > » (n being interpreted so ^ na ^ 
as time) , and (B) is the well-known Liebmann method. 

We assume now that (3) is replaced by the stronger 
condition 

F z >v>0 (49) 



m n <p n ~ 1 m 1 , 



and prove that procedure (B) converges. The proof 
for (A) is practically the same. 



and for #>1, 

max fi j0 (w+ ^-0 (w) |< §J m v < ^ 

v =n+l 1 — P 



By (48) 



which implies the existence of 

0=lim <t> {n) . 



^+»-^>=-g (Ipy^^-??,-^, (50) Also, by (46), (48), 



where the asterisk indicates that the partial deriva- 



# +, -*=s(fp) (3W-**), 



(51) 



(52) 
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where the double asterisk indicates that the partial and satisfying the following conditions: 
derivatives are evaluated at a point (f 1 , . . ., f 4 ), 



A/,[0]=A[<£] if <t>{x,y) is a polynomial of second 

degree, (59) 

p"(x,y)>0for(r,«)^(0,0), (60) 

o<^-< S P ff (*,y)£*i. (6i) 

Applying (59) to = 1 we have that 

P 00 (x,y)=- S P TS (x,y)<0. (62) 

(r,s) = (0,0) 



r=T«W+(l-T0*lr, 0<T,<1. 

Hence, setting 

/i»+i=maxo> ( » +1 >-0|, (53) 

we have that /x n+ i <pix n} so that 

Hn<P n n>. (54) 

Since /z can be estimated by means of (31), we can 

determine a priori how many steps are needed in 

order to obtain <t> with a desired degree of accuracy. ? +<?$ has continuous derivatives of the second 

Now we drop assumption (49) and set (for „>0) ™jer whose moduli of continuity do not exceed 

LftoJ = £„[<*>(]+#*• 

Let Let <j>^ be the solution of the boundary value 
problem 

£H4>J = 0, i=l, . . ., N; 

4>t=Vi, J=N+1, . . ., M. 



|A A M-AM|<2^ 1 o,< 2 >(V2^), 



(63) 



and if 4>(x,y) has partial derivatives of the third 
order bounded in modulus by m (3) , then 



lAM-AMl^fK^h. 



(630 



This solution exists by virtue of the preceding result. Under the hypotheses that the fourth order deriva- 
Applying (33) with L h replaced by Lf\ we conclude tives exist and are bounded by m (4) in modulus this 
that |0 ( ^ | <K, so that \L h U^]\ = \rj^ v) | <Krj. By (31) can be improved to 



so that 



max OA |0 ( «'> -^"'l^CKW+rj"), 



</>=lim ^ 



exists. It is clear that ^ is a solution of P h 
by (31), 

maxoj^ — cj>\<CKr]. 



{AM-AMlKp^mMh*, (63") 

provided A,, satisfies the additional condition 

A1 A h [<j>] =0 if <Kx 9 y) is a homogeneous 

iUS0 > polynomial of degree 3. (64) 

^ ' All this follows easily from Taylor's theorem. 



Since <£ (l,) can be computed by procedures (A) or (B), 
we have an effective method of solving P h in all cases. 

7. Extension 

Let a(x,y), b(x,y), c(x,y) be continuous functions 
defined in Q. 

ac-6 2 >0. (56) 

Set 

A[cj)]=a<t> xx +2b<l) X y J rC(l> yv , (57) 

Lfo] - Afo] - F(3,y,*,*tA) • (58) 

With this new operator L we can state problem P. 



Once A>> has been chosen, it is easy to find differ- 
ence operators 

1 i 

Dd,z[<t>]=T S a rs (x,y)<t>(x+rh,y+sh), 

to r,s=-t 



Dh, y [<t>]=T S T rs (x,y)<t>(x+rh,y + sh), 

to r,s = -t 

satisfying the conditions 

„«>(*,*) -T*»(*,y)«0, 
\* T, (x,V)\y(x,V)\<K2\p"(*,V)\, 



Motzkin and Wasow* showed that (t being a 
sufficiently large positive integer depending on A) , , . ,, ,, 

A can be approximated by a difference operator and approximating ^ <$> y in the sense that 

A fc (fc>0) defined by 



(65) 
(66) 



1 



AftfoH-o S P TS (x,y)<t>(x+rh,y+sh), 

to T,S=-t 



•J. Math. Phys. 31 (1953). 



Dn x[4] — 4>x D h y [<t>] = <j) y for every linear function <t>(x,y). 

(67) 
We have that 

IDnM-^DnM-^l^tK^^i^th) (68) 
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if <t>(x,y) has continuous partial derivatives with 
moduli of continuity bounded by co (1) (77) , 

\D h , x [<t>]-4>x\,\DnM-^\ <2t*K x K 2 m,v>h (680 

if the second derivatives exist and are bounded by 

m (1) in absolute value. 

If 

A*,rr[<£]=A*,j0]=O for every homogeneous quadratic 
polynomial 4>(x,y), (69) 

then, under the hypothesis of (63'), 

\D h A4>]-<t>^\D h M-<i>y\ <^K x K 2 m^h\ (68") 
We approximate L by the difference operator 

L h [4>] = AnW ~ F{ x,y,<t>,D h >x fo] ,D h tV [<fi } . (70) 

Let P be the point (x ,y ). Weorder the T=4tt(t+1) 
points lx +rh,y +sh) y where r,s = 0,±l, . . ., ±t, 
and (r,s) 9^ (0,0) lexicographically with respect to 
(r,s) and denote them by P i, . . ., Pot- These 
points shall be called the ^-neighbors of P . The 
new definition of neighbors leads to a new definition 
of the boundary £l' h of a lattice domain Q h , and having 
this new definition, we can state the boundary value 
problem P h for the operator L h defined by (70). 

We make now two remarks concerning the linear 
difference operator 

W*<] = A*[0i] + a*D».x[0j + j8|D*,ir[*J-T^* (71) 



so that l h can be written in the form (73) with 

\ ip =p^ 8) '(x i9} y iv )— h<r ir > s) *(x iv ,y iv )— hr {r ' s) >>(x iv ,y iv ), 

where P iv ={x i -{-rh,y i J rsh), (r,s) = (r,s) v . By virtue 
of (61), (66), condition (72) implies (74). 

(2) Let l h denote the operator (71) with 7*^0. 
There exists a continuous function $(x,y) defined in 
12 and positive constants h 2 , C such that for Q<^h<Ch2 



subject to conditions (13), (14). 
(1) Set 

hl== 2K 2 
Then for 0<A<Ai 



with 



k[*i]=jp {2jx ,r *ir-*«Z}x , 'J-7«** 



\'">0, Z> f *> - 



(72) 

(73) 
(74) 



nmx- Q $(x,y) <c 
mmzT h [$(x,y)\ 



(75) 



In fact, assume (without loss of generality) that 
ft is contained in the domain 0<x<d, and set 
$ = ^ M >0. By (63 r ), (680 we have that 

UHx,y)]>a{x,y)^ x -^K x ^ d h 

—A i ie iiX -2A2t 2 K l K 2 ii 2 e iid h . 

If we first choose a fixed ju such that 

M 2 mino a (»,i/) — /xA>l, 

we can determine a h 2 such that for h<ih 2 ,l h [$>(x,y)]<^%. 
Then (75) holds with C=2e» d . 

Using these two remarks and assuming that 
/K^min (hi,h 2 ), all considerations of sections 3 to 6 
can be repeated with only minor and obvious modi- 
fications. Inequality (16) holds with the C deter- 
mined above, and so does (31). This inequality 
implies the uniqueness and existence of a solution of 
P h . This solution can be computed effectively by 
the method of section 5, and its deviation from the 
solution of P can be estimated as in section 5, using 
inequalities (63), (68). 

We note that the case in which eq (1) is approx- 
imated not by (2) but by a more sophisticated 
difference equation is included in the preceding 
discussion. 



In fact, applying (57) to <£=1, we see that 
2cr-(^)^2r-(^)^0, 



Los Angeles, October 9, 1952. 
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